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We present the results of the tree-level calculation of inclusive two-hadron production in 
electron-positron annihilation via one photon up to subleading order in 1/Q. We consider the 
situation where the two hadrons belong to different, back-to-back jets. We include polarization 
of the produced hadrons and discuss azimuthal dependences of asymmetries. New asymmetries 
are found, in particular there is a leading cos(2</>) asymmetry, which is even present when 
hadron polarization is absent, since it arises solely due to the intrinsic transverse momenta of 
£S| ■ the quarks. 
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Three basic hard scattering processes in which the structure of hadrons is studied with electroweak probes 

are (semi-)inclusive lepton-hadron scattering, the Drell-Yan process and inclusive hadron production in e + e~ 

annihilation. In this paper we will focus on the latter, where we restrict ourselves to photon exchange. The 

^ i| (timelike) photon momentum q sets the scale Q, where Q 2 = q 2 , which is much larger than characteristic 

hadronic scales. 

The inclusive lepton-hadron scattering, generally known as deep inelastic scattering (DIS), is the most studied 
process from the theoretical as well as from the experimental side. From the theoretical point of view DIS can 
be described by using the operator product expansion (OPE), within the context of Quantum Chromodynamics 
(QCD). This allows to relate moments of structure functions to (Fourier transforms of) hadronic matrix elements 
C^ ■ of local operators. In the (QCD improved) parton model, i.e., to leading order in 1/Q, the structure functions 
can be expressed as sums of so-called distribution functions (DF's). The OPE thus gives information on moments 
of the DF's. However, one can also describe a scattering process directly in terms of the DF's themselves, which 
are (Fourier transforms of) hadronic matrix elements of non-local operators jl|,|2|] . For treating the non-leading 
orders in 1/Q, Ellis, Furmanski, Petronzio (EFP) |1 have developed a formalism for unpolarized DIS, using 
non-local operators, following ideas by Politzer 8. The extension to polarized DIS was done by Efremov and 
Teryaev S. At tree level, i.e., order (a s )°, EFP have shown the equivalence of their formalism to the OPE 
approach for (1/Q)™ power corrections. 

The non-local operators consist of quark and gluon fields. The quark (gluon) DF's are functions of the light- 
cone momentum fraction x = p + / P + of a quark (gluon) with momentum p in a hadron with momentum P. 
For multi-parton DF's, which show up beyond leading order, the functions depend on momentum fractions of 
several partons [H^]. Often the number of partons can be reduced with help of the QCD equations of motion. 

The EFP approach was extended to other processes in analogy to the DIS description, in particular, to those 
ones where the OPE is not applicable (for an alternative approach based on a non-local operator expansion 
see 0). The Drell-Yan process has been studied to leading order by Ralston and Soper || and partially to 
subleading order in |p|dTl|. Here the cross-section is a sum of products of two DF's. Recently, the complete 
tree-level result up to order 1/Q for semi- inclusive polarized lepton-hadron scattering was published |l2|]. In 
this case one needs to include so-called fragmentation functions (FF's) p~^,0] and the cross-section is a sum of 
products of a DF and a FF. 



The present paper focusses on the third process of interest: inclusive two-hadron production in e + e~ annihi- 
lation up to order 1/Q, where the two hadrons belong to different, back-to-back jets. The cross-section involves 
products of FF's, the number of which is larger than the number of DF's, due to the appearance of so-called 
time-reversal odd structures, which arise since time reversal does not give constraints in this case. Interesting 
features like the Collins effect E3] show up, as do other new asymmetries. 

The DF's and FF's are essentially non-perturbative objects and must be determined by experiment or cal- 
culated, for instance with the help of models. The idea is that these functions are universal and once they are 
measured in one process they can lead to predictions in others. The DF's and FF's parametrize the structure of 
a hadron, e.g. the spin structure, and the leading order DF's and FF's have simple probabilistic interpretations. 
With the help of symmetry properties one can first of all determine the set of DF's and FF's in which a process 
can be expressed. By giving cross-sections and asymmetries in terms of these functions one can deduce how to 
measure and separate them. 

From the experimental side most well-known are the leading order unpolarized and polarized DF's called /i 
and <7i, respectively, Q and the FF D\ (l(|. For subleading order the DF g 2 has been measured ]l7] , but 
still with rather large errors. The experimental knowledge on FF's is much smaller than that on DF's. We 
obtain the most general expression for the cross-section in terms of as yet unknown FF's, in order to find out 
which other functions might be experimentally accessible in the near future and where to look for them. Some 
of the unknown functions have been modelled p8[ , which results can be used to estimate the magnitudes of 
asymmetries from our expressions. 

Some asymmetries, like the Collins effect, are leading effects, not suppressed by powers of 1/Q, but one needs 
to include intrinsic transverse momentum in the DF's and FF's p]Mp|. Intrinsic transverse momentum plays a 
crucial role in two-hadron processes involving two soft, non-perturbative, parts, since the transverse momenta 
are linked by momentum conservation. Those DF's or FF's which, as functions of transverse momentum, 
would not contribute in one-hadron processes, will show up in these two-hadron processes. The idea that 
intrinsic transverse momentum always gives rise to suppression is incorrect, although it is the case in one-hadron 
processes. In the pioneering work |19|-|2l| on azimuthal dependences due to intrinsic transverse momentum, all 
effects arc found to be suppressed by at least 1/Q. 

On the formal theoretical side there remains the open issue of how to proof factorization for the process under 
consideration. For the case of back-to-back jets there exists a proof, but no higher twist effects are included 
p2| . We do not expect polarization to be a problem for factorization (see [£3| ) . For the case of the Drell-Yan 
process, which is very similar to the process under consideration, arguments have been given why factorization 
holds for the first non- leading power corrections |2jj2^], which is consistent with its known failure at order 1/Q 4 
p(i| p8[ . We expect factorization to hold also for the case at hand. 

We will not concern ourselves with these problems here, even though factorization is needed to ensure uni- 
versality of the FF's, and restrict the discussion to tree-level. It represents the extension of the naive parton 
model to subleading order and shows the dominant structures to be expected in the cross-section, although 
QCD corrections, such as Sudakov effects p9| , may affect the magnitude of the asymmetries. At tree-level the 
only QCD input at order 1/Q (apart from the Feynman rules) is the use of the equations of motion which 
ensures the electromagnetic gauge invariance. 

The outline of this paper is as follows. In Section 2 we present the formalism of the e + e~ annihilation 
process, with emphasis on the kinematics. Section 3 contains the analysis of the soft parts of the process, in 
particular the fragmentation functions are studied. This is followed by the details of the complete tree-level 
calculation of the hadron tensor up to subleading order in Section 4, the result of which is given in Appendix 
B. In the three following sections we investigate special cases, which give more insight than the full result and 
are useful from a practical point of view. In Section 5 we discuss the result after integration over the transverse 
momentum of the photon. In Section 6 we study the differential cross-section, i.e., not integrated over transverse 
photon momentum, but restricted to leading order and the case where only one hadron is polarized (the case 
of two polarized hadrons is given in Appendix C). In Section 7 this is compared with the integrated cross- 
section weighted with factors of the transverse momentum of the photon. Finally, the results are summarized 
in Section 8. 



II. KINEMATICS 

We consider e~ + e + — v hadrons, where the two leptons with momenta I and V annihilate into a photon with 
momentum q = I + V , which is timelike with q 2 = Q 2 — > oo. Denoting the momentum of outgoing hadrons by 
Ph (h = 1) 2, . . . ) we use invariants z/j = 2Ph ■ q/Q 2 . The momenta can also be considered as jet momenta. 



We will consider the general case of polarized leptons with helicities ±A e and production of hadrons of which 
the spin states are characterized by a spin vector Sh (h = 1, 2, . . . ), satisfying S^ = — 1 and Ph • Sh = 0. In 
this way we can treat the case of unpolarized final states or final state hadrons with spin-0 and spin-1/2. We 
will work in the limit where Q 2 and Ph ■ q are large, keeping the ratios Zh finite. 

The square of the amplitude can be split into a purely leptonic and a purely hadronic part, 

\M\ 2 = ^L^H^, (1) 

with the helicity-conscrving lepton tensor (neglecting the lcpton masses) given by 

L^(l, I'; A e ) = 21X + 2U; - Q 2 g^ + 2i\ e e^paM" '. (2) 

For the case of two observed hadrons in the final state, the product of hadronic current matrix elements is 
written as 

H^(P X ;P 1 S 1 ;P 2 S 2 ) = (0\J^(0)\P x ; P^; P 2 S 2 )(P X ; P 1 S 1 ; P 2 S 2 \JA0)\0), (3) 

where a summation over spins of the unobserved out-state is understood. The cross-section (including a factor 
1/2 from averaging over incoming polarizations) is given by: for 2-particle inclusive e + e~ annihilation 

<PPi<PP 2 4Q« ""^ ' W 

with 

1 I cfiP 
W^(q; Pi Si; P 2 S 2 ) = -— < J (2 ^3 2 p (2^) 4 * 4 (g -Px-Pi- P2)H^(P x ;P l S l ;P 2 S 2 ), (5) 

for 1-particle inclusive e + e~ annihilation 
with 

i r d 3 p x 4 . 



W^(q;P h S h ) = — - \ " {2vf5\q-P x -Ph)(0\MO)\Px;PhS h )(Px;PhS h \J v (0)\0), (7) 

(Ztt) J {ZTr)- i 2P x 

and for the totally inclusive annihilation cross-section the well-known result 

4tt 2 a 2 
o~(e e — > hadrons) — — — — — L llll R >lv , (8) 



with the tensor R^ u given by 

R,Aq) = j j0^(^) 4 s 4 (q- Px)(o\MO)\p x )(p x \MO)\o) 

= [ d*x e* x {OW^xlMOWO). (9) 



Recall that the totally inclusive cross-section is directly related to the vacuum polarization. Also note that 
the totally inclusive process is short-distance dominated, whereas the 1-particlc inclusive case is light-cone 
dominated, but only in the former the OPE can be applied. 

In order to expand the lepton and hadron tensors in terms of independent Lorentz structures, it is convenient 
to work with vectors orthogonal to q. A normalized timelike vector is defined by q and a normalized spacelike 
vector is defined by P M = P^ — (P ■ q/q 2 ) q^ for one of the outgoing momenta, say P 2 , 

a'" 
*»- W P% = 2^--%-. (11) 
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FIG. 1. Kinematics of the annihilation process in the lepton center of mass frame for a back-to-back jet situation. 
Pi is the momentum of a jet or of a fast hadron in a jet, Pi is the momentum of a hadron belonging to the other jet. 



Note that we have neglected 1/Q 2 corrections, as we will do throughout the paper. Such corrections arise 
among others from the hadron masses Mh, so-called target mass corrections or kinematic power corrections. 
Vectors orthogonal to z and i are obtained with help of the tensors 



gtu> = gi iu _ pq, + 2»z u , 



pv _ iivpaj. - 

r - J_ c L p^(T 



(Pi ■ q) 



^pvpo p „ 



(12) 
(13) 



For instance, using the other hadronic momentum Pi , one obtains P{Y = g^ P\ v (see Fig. [l]). We define the 
normalized vector h^ = P[Y/|Pi±| and the second orthogonal direction is given by e^h v . We use boldface 
vectors to denote the two-dimensional Euclidean part of a four- vector, such that Pu_ • Pi_l = — Pi± ■ Pi±- In 
the calculation of the hadron tensor it will be convenient to define lightlike directions using the hadronic (or 
jet) momenta. Consider two hadronic momenta Pi and P2 not belonging to one jet (i.e., their dot product 
Pi • P 2 is of order Q 2 ). The momenta can then be parametrized using dimensionless lightlike vectors n + and 
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and n + 
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(14) 
(15) 



^vr- + vT+ + 9 - 



(16) 



where Q 2 = 



Qip with q 2 -, = —Q 2 ^. We will use the notation p — p ■ n + and p + = p ■ n_ for a generic 



momentum p. For the case of two back-to-back jets Q\ <§C Q 2 and up to Q^/Q 2 , which we neglect, one has 
Q = Q, (1 = zi and C2 = z 2 . If momentum P 2 is used to define the vector z M , then 



pi' 

Ml 



-ziq£. 



(17) 



Vectors transverse to n + and n_ one obtains using the tensors 



9t 



^pv 



{p v\ 



pl^vpcr 



(18) 
(19) 



where the brackets around the indices indicate symmetrization. Note that these transverse tensors are not 
identical to the perpendicular ones defined above if the transverse momentum of the outgoing hadron does not 
vanish. The lightlikc directions, however, can easily be expressed in t, z and a perpendicular vector, 
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(20) 
(21) 



showing that the differences are of order 1/Q. Especially for the treatment of azimuthal asymmetries, it is 
important to keep track of these differences. 

In summary, we use two sets of basis vectors, the first set constructed from the photon momentum (q) and 
one of the hadron momenta (P2), the second set from the two hadron momenta (Pi and P2). The respective 
frames where the momenta q and P%, or P\ and P2, are collinear are the natural ones connected to these two 
sets. In the first Pi has a perpendicular component P±±, in the second q has a transverse component q T . One 
can get from one frame to the other via a Lorentz transformation that leaves the minus components unchanged 



III. CORRELATION FUNCTIONS 

In this section we discuss the relevant 'soft' hadronic matrix elements that appear in the diagrammatic 
expansion of a hard scattering amplitude. Assuming the two hadrons to belong to two different jets we encounter 
two types of soft parts in the process under consideration: one describes the fragmentation of a quark into a 
hadron plus a remainder which is not detected and the other describes the similar fragmentation for an antiquark. 
Up to order 1/Q the quark fragmentation is described with help of two types of correlation functions: the quark- 
quark correlation function A(Pi, S±;k) B and the quark-gluon-quark correlation function A^(Pi, Si; k, k\) 
(Fig- 1): 
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FIG. 2. Correlation functions A and A a- 



A u -(P,Si;fc)=^^ i y'd 4 x e * a; (0|^(x)|Pi,Si;A)(P,S 1 ;A|^(0)|0), (22) 

A^ u .(P,Si;A:,fci) = ^^y'd 4 xAe Jfc - y+lfcl - ( ^ y ^O|^(a ; ) 5 ^(y)|Pi,S 1 ;A)(P,Si;X|^.(0)|0), (23) 

where k 7 k\ are the quark momenta and an averaging over color indices is understood. If one chooses the gauge 
A~ = only a transverse gluon is relevant. In fact, in a calculation up to subleading order, we only encounter the 
partly integrated correlation functions J dk + A(Pi, Si; k) and J dk + d A k\ AJ(Pi, Si; k, fci), which are functions 
of k~ and k T only. Note that the definition of A^ includes one power of the strong coupling constant g. 

The above matrix elements as functions of invariants are assumed to vanish sufficiently fast above a char- 
acteristic hadronic scale, which is much smaller than Q 2 . This means that in the above matrix elements 
k 2 , k ■ P\ ^C Q 2 . Hence, we make the following Sudakov decomposition for the quark momentum k: 
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Similarly, we decompose the spin vector Si: 

\\Z\Q 
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(24) 



(25) 



with for a pure state A^ + S\ T = 1. In the approximations the + components (oc 1/Q) are neglected, as these 
are irrelevant compared to the + components of the momenta in the hard part (oc Q). 

The Dirac structure of the quark-quark correlation function can be expanded in a number of amplitudes, i.e., 
functions of invariants built up from the quark and hadron momenta, constrained by hcrmiticity and parity 
P,12[. Here we directly integrate the correlation function over fc + , which up to order 1/Q can be parametrized 
as follows: 
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(26) 



where the shorthand notation G\ s stands for the combination 



G u (z, k T ) = Ai Gil + G 1T 



(fc, 



'IT) 



Mi 



(27) 



etc. The functions E, Di, ... in Eq. ( |26|) and Gil, Git, ■ ■ ■ in Gi s , . . . are fragmentation functions. One wants 
to express the fragmentation functions in terms of the hadron momentum, hence, the arguments of the frag- 
mentation functions are chosen to be the lightcone (momentum) fraction z = P^ /k~ of the produced hadron 
with respect to the fragmenting quark and k' T = —zk T , which is the transverse momentum of the hadron in 
a frame where the quark has no transverse momentum. In order to switch from quark to hadron transverse 
momentum a Lorentz transformation leaving k~ and P^ unchanged needs to be performed. The fragmentation 
functions are real and in fact, depend on z and k' T 2 only. 

Inverting the above expression, the fragmentation functions appear in specific Dirac projections of the corre- 
lation functions, integrated over k + : 



A^( Zl k T ) = — f dk+ Tr(AT) 
Az J 



=Pj /z, k T 



E 

X 



dx + d 2 x 



T ik-x 



4z (2tt) s 



■ x Tr(0\^(x)\P 1 ,S 1 ;X)(P 1 ,S 1 ;X\^(0)T\0) 



x~=0 



for which we can distinguish the leading fragmentation functions: 



(28) 



Ah-! (z, k T ) = Dl (z, k' T 2 ) + e l^lll dj. Axi fc -), 

A^'^(z,k T ) = G ls (z,k T ), 

A^-t*] (z, k T ) = S\ T H lT (z, k' T 2 ) + || Ht.{z t k T ) H 



f J k 
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furthermore we obtain sublcading projections (i,j are transverse indices): 

Mi 



A^(z,k T ) = -^E(z,k' T 2 ), 

A^(z,k T ) = ^DHz,k' T 2 ) 
A^](z,fc T ) = ^£ s (z,fc T ), 



^|^i(z,fc T 2 ) 



^_ 1Tj D T (z,k' T 2 ) 



(29) 
(30) 
(31) 

(32) 

(33) 
(34) 



A h^} (2j kr) = ^it G ^ ^ + h_ G ± (Zi Mj (35) 



p- ^ , , p. 

Al-^1 (z, k T ) = S lTk° T -k l T S{ T H ^ fe ^ 2) + M0_ H{ ^ ^ (36) 

A^ + ^(z,fc T ) = ^iJ s (z,fe T ). (37) 

We identified leading and subleading functions, which in principle start contributing at order 1 and l/Q, 
respectively. The order at which a function first can contribute depends on the power of M\/P^~ in front of the 
function as it appears in the projections. Each factor M\j P^~ leads to a suppression with a power of M\/Q in 
cross-sections. We will refer to the function multiplying a power (Mi/Pr~)*~ 2 as being of 'twist' t. We note 
that this notion of twist, in analogy to the k T integrated case M], is related but not equal to the one used for 
local operators in the OPE. 

The naming scheme is as follows. All functions obtained after tracing with a scalar (1) or pseudoscalar (175) 
Dirac matrix are given the name £.., those traced with a vector matrix (7^) are given the name D , those 
traced with an axial vector matrix (7^75) are given the name G.. and, finally, those traced with the second rank 
tensor {io^^s) are given the name H . A subscript 1 is given to the leading functions, subscripts L or T refer 
to the connection with the hadron spin being longitudinal or transverse and a superscript _L signals the explicit 
presence of transverse momenta with a non-contracted index. In the literature sometimes the fragmentation 
functions are denoted by lower-case names, but supplemented by a hat (e,g,h), with the one exception that 
D is named /. We note that after integration over k T several functions disappear. In the case of A^ 75 ' 
and A^ 75 ! specific combinations remain, namely H\ = Hit + {k T /2Mf) H^ T and Gt = G' T + (k T /2Mf) G^, 
respectively. 

The choice of factors in the definition of fragmentation functions is such that J dz d 2 k' T D\ (z, k' T ) = Nh, where 
Nh is the number of produced hadrons. The twist-two fragmentation functions have natural interpretations 
as decay functions. The projection A^ 7 1 is (after proper normalizing) the probability of a quark to produce 
a spin-1/2 hadron in a specific spin state, A^ 7 75 1 is the difference of the probabilities for a chirally right and 
chirally left quark to produce such a hadron, while A^ 75 ' is the difference of opposite transverse spin states 
(along direction i) of a quark to produce such a hadron. 

Note that the decay probability for an unpolarized quark with non-zero transverse momentum can lead to 
a transverse polarization in the production of spin-1/2 particles. This polarization is orthogonal to the quark 
transverse momentum and the probability is given by the function D^ T . In the same way, oppositely transversely 
polarized quarks with non-zero transverse momentum can produce unpolarized hadrons or spinless particles, 
with different probabilities. This difference is described by the function H^-, which is the one appearing in the 
so-called Collins effect E3|, which shows up as a single transverse spin asymmetry in semi- inclusive DIS, and 
arises due to intrinsic transverse momentum. 

The functions D^ T and H^ are examples of what are generally called 'time-reversal odd' functions. This 
somewhat misleading terminology refers to the behavior of the functions under the so-called naive time-reversal 
operation T/v fl32] , which acts as follows on the correlation functions: 

A(P 1 ,S l ;k) 2U ( l5 CA(P 1 ,S 1 ;k)C^ 5 y (38) 

where k — (fc°, — k), etc. If T/v invariance would apply, the functions D^ T , H^~, Dj;, Dt, El, Et and H would 
be purely imaginary. On the other hand, hermiticity requires the functions to be real, so these functions should 
then vanish. 

The operation T/v differs from the actual time-reversal operation T in that the former does not transforms 
in into owi-states and vice versa. Due to final state interactions, the out-state \P±, S\; X) in A(Pi, S\) k) is not 
a plane wave state and thus, is not simply related to an in-state. Therefore, one has T/v 7^ T and since T itself 
does not pose any constraints on the functions, they need not vanish. 

In the analogous case of distribution functions, which are derived from matrix elements with plane wave 
states, T = T/v and therefore there are no 'time-reversal odd' distribution functions. 

The H and E functions (but not E s ) are called chiral-odd functions because they are non-diagonal in the 
chirality basis, so they arise either accompanied by a quark mass term or by another chiral-odd function, such 
that the product is again chiral-even [p3p|] . 

The quark-gluon-quark correlation functions can be expressed in terms of the quark-quark correlation func- 
tions with help of the classical equations of motion (e.o.m.). These can be used inside hadronic matrix elements 
M. If we again define Dirac projections: 



A^ [r] (z,fc T )= A. /'rffc+d 4 fc 1 Tr(A^r) 

= £ 



4z (2tt) 
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we find as a consequence of the e.o.m.: 
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where the functions indicated with a tilde (7J, E, . . . ) differ from the corresponding twist-3 functions (H , E, 

. . .) by a twist-2 part, namely 
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(44) 
(45) 

(46) 

(47) 
(48) 

(49) 

(50) 
(51) 
(52) 

(53) 



We have included Gt in this list since it is relevant for fe^-integrated functions and note that Gt = G' T + 
(k T /2Al 2 ) Gji. The functions in Eqs. (40) to (42) are interaction-dependent and vanish for the case of a quark 



fragmenting in a quark (as can be checked with the help of Appendix A). Note that the time-reversal odd twist- 
2 functional D^ T and H^ are in fact interaction-dependent. Their presence is due to final state interactions 
of the produced hadrons, which after all are strong interactions. The separation of twist-3 functions in this 
way is analogous to the case of the distribution function qt = gx + .92 , and the twist-2 parts could be called 
Wandzura-Wilczek parts [p4[. 



x The arbitrariness in the definition of Dt and H in Eqs. (k6|) and fed) is fixed by the requirement that the functions 
D± T and H^ do not appear in the integrated versions of Eqs. mu) to (B4) . 



For the fragmentation of an antiquark most things are analogous to the quark fragmentation. The major 
difference in our case is that the role of the + and — direction is reversed. We will denote the antiquark 
correlation functions by A(P2, S 2 ;p). These should be defined consistently with the replacement ip — > ip c = Ctp T , 
or A I r l = A c I r l for T = 7^, icr^^s, 275 and A ^ = — A c I r l for T = 1, 7^,75, where we have defined the projections 
as: 



A [rl (z,p T )= -^jdp-Tr(Ar) 



(54) 

P+=P+/Z, P T 



where we make the following Sudakov decomposition for the antiquark momentum p: 

P=— =n+ + n_+p T » -P 2 +p T . (55) 

zV2 z 2 Q\/2 z 



Similarly, we decompose the spin vector S\ 



2- 



X 2 z 2 Q X 2 M 2 A 2 . . 

S2= A^^^-^Q^"- +S2r "^ P2 + ^' (56) 

with for a pure state A 2 + S\ T — 1. The antiquark fragmentation functions are denoted by D 1 (z, p^ 2 ), . . ., with 
p' T = —zp T , in full analogy to the quark fragmentation functions. The antiquark fragmentation functions are 
obtained from 

S y (P 2 ,5 2 ;p) = St^u /<*** e^ x (0|? J (0)|P 2 , > S 2 ;X)(P 2 ,5 2 ;X|^(x)|0). (57) 

Although the antiquark-gluon-antiquark correlation functions are straightforwardly defined, we still give here 
the relations which follow from the e.o.m., since these differ non-trivially from those for the quark-gluon-quark 
correlation functions. We will not use tilde functions here, since in the non-symmetric frame in which we will 
express the hadron tensor (non-symmetric between quark and antiquark fragmentation part), they do not show 
up in a natural way. 

A A a = £ t a a = *(— E-mD 1 + 1 — E-i—^ E 1 

1^=1. . m — _i_ 



e%p Ti S 2Tj I H 1T - -E T +i — D 1T 1 , (58) 

A , ;; ; = -^E s + mG ls + i ^E s + (p T ■ S 2T )Ext + %- K, (59) 

z z M 2 

A A -ze T M A A/3 =-p T \-D - D 1 -i — E l I - —^p Tl S 2T] D 1T 

■ «/3 ( l n^ n m Tr 1 , ■ X2 n ± 

- k t PTfi I ^ G s ~ Gls ~Jf H ^ + l Y L 

— te 7 rb 2Tf3 I — ^- G T — uiHxt + i ~r- l>t I • (60) 

Until now we have not commented on color gauge invariance of the correlation functions. As given above 
they are gauge-invariant quantities displayed in a specific gauge. In general, one has to include path-ordered 
exponentials, in order to compensate for the gauge non-invariance due to the non-locality of the operators. Such 
a link operator is of the form: 



C(0,x)=Vexp(-ig dz* A^ 



(z) . (61) 



At this point we assume that matrix elements with multiple A -gluon fields in A^ , A AA , . . . (multiple A + 
gluon fields in A^, A^, . . .) will combine into an appropriate link operator with path along the + direction (- 



direction) in A (A) (cf. [jlCj). For the fc T -dependent functions which involve transverse separations, the path 
from the point to x in A will run along the + direction via x + — oo. The transverse part of the path, which 
is at co, does not contribute, since matrix elements are assumed to vanish there. 

There remains one issue to be addressed, namely the explicit A T in A^ is not gauge invariant. Nevertheless, 
using the covariant derivative, one can express this A T in terms of D T and <9 T , such that: 



Af\z,k T ) 



v«PT 



z,k T )-k a A [r] {z,k T ), 



where 



Y — 



d^x e l "- x (0\ij i (x)id^\P 1 ,Sv,X)(P l ,Sv,X\ilj j (0)\0) = k' i Ai^PuS^k). 



(62) 



(63) 



Hence, we see that (again after inclusion of a link operator, assumed to arise from A^, A^TJ^, . . .) A\ is a color 
gauge invariant quantity and therefore, so are the interaction-dependent parts of the twist-three fragmentation 
functions. 



IV. THE COMPLETE TREE-LEVEL CALCULATION 

Up to order 1/Q there are five tree-level diagrams to consider. The simplest diagram (Fig. |3|) involving only 
quarks contributes at order 1 and 1/Q, the other four (Fig. 0) involve one gluon which connects to one of the 




2 2 

FIG. 3. Quark diagram contributing to e + e~ annihilation in leading order. There is a similar diagram with reversed 
fermion flow. 



two soft parts. Note that one power of the coupling constant is included in the definition of the soft part, such 
that the diagrams are of order (a s )°. 

The momentum conserving delta-function at the photon vertex is written as (neglecting 1/Q 2 contributions) 

8\q -k-p) = 6(q + - p+) 5(q- - k~) 8 2 {p T + k T - q T ), (64) 

fixing P 2 + /z = p + = q + = P 2 + '/ ' z<x and P^ / z = k~ = q~ = Pj - / ' z\. Eq. ( p4[) shows why only the k + and 
p~ -integrated correlation functions are relevant. Note that the quark transverse momentum integrations are 
linked. The five diagrams lead to the following expression for the full result up to order 1/Q: 



W 



--^/♦-^^^^-^(StoyW) 



- Tr ^A a A (p)rA(k) la -^±=rj Tr (J 7 „A^(p) 7o ) ^±±= la ^k)^ 
+ Tr (A(p)r ( 7 oA^(fc) 7o ) 7^^7 Q ) + Tr (A{p) la ^=^ A£(fc)7 v ) 



p+k- 



(65) 
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^p, 



kJX/X, 

p, L 




,lk 



pi- 




lf\J\, 



p|- 




^•' '^ 



"|p 



_q_ k| 

p|- 




r \rvr 



2 2 r 2 2 

FIG. 4. Diagrams contributing to e + e~ annihilation at order 1/Q. 

The factor 3 originates from the color summation. We have omitted the flavor indices and summation; further- 
more, there is a contribution from diagrams with reversed fermion flow, which results from the above expression 
by replacing fj, <-> v and q — > —q. 

In the expression the terms with fi± arise from the fermion propagators in the hard part neglecting contri- 
butions that will appear suppressed by powers of Q 2 , 

4 - j>\ + m ^ (g+ -p\ + Yr _ 1_ _ ^+ 
(q - pi) 2 - to 2 ~ 2( g + - pi+)g- ~ 2q~ ~ Q^ 

#i - j+m (fci~ -<?~h + _ 7 + = J- , fifi ^ 

(fc!- g ) 2 -m 2 ~ -2{k 1 --q-)q+ -2q+ Q^ K ' 

where the approximate sign holds true only when the propagators are embedded in the diagrams. The quan- 
tity A^(k) arises from integrating out the second argument of A^(fc,/ci) instead of the first which yields the 
combination 7oA^ (fc)7o: 

j ^ A Stf( P i' S i; fc '*i) = E(2^ j ^ d 4 x e ik - x (0\Mx)9A^(x)\P 1 ,S 1 ;X)(P u S 1 ;X\iJ j (0)\0) 

= A% ij (P 1 ,S 1 ;k), (67) 



= ( 10 Af 10 ) ij (P 1 ,S 1 ;k) 

and similarly for A A (p) and 70 A ^ (p)jo- To deal with these combinations one can use the relation: 

,[r] 



(7oA^ f 7o 



A° 



(68) 



(69) 



and a similar one for 70 A^ (p)7o- 
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To obtain the expressions for the symmetric and antisymmetric parts of the hadron tensor we expand all 
vectors in A, A, AJ4 and A^ in the perpendicular basis (t, z and _L directions). In particular, we reexpress the 
transverse vectors k T , p T , S 1T and S 2T in terms of their perpendicular parts and a part along t and z. For this 
we need 



IIV 

9t 



up V 

9T9t p 



1*4 



(70) 



We will refer to the perpendicular projections as k±, etc. (instead of the fully logical name, which would be 



k T± ). Thus 



/Cj_ g^_ KTj, Krp 



V T«/ 






(*" + «"), 



(71) 



and similarly for pj_, S\j_ and S2j_. We note that for these four vectors with this definition the two-component 
perpendicular parts are the same as the two-component transverse parts, i.e., k± = k T , Si± — S 1T , etc. 

The full expressions for the symmetric and antisymmetric parts of the hadron tensor (expressed in the 
perpendicular frame defined in section 2) are given in Appendix B. We note that the expressions are not 
symmetric in the interchange of the hadrons 1 and 2, because the choice of perpendicular direction (P21. = 0) 
is non-symmetric. 

The cross-sections are obtained from the hadron tensor after contraction with the lepton tensor 



L? v = Q 2 



- (1 - 2y + 2y 2 ) ft + 4y(l - y)z»z" 

-4tf(l - y) {ln v + \ <?r) - 2(1 - 2y)^y~^y) z { H^ 

+i\ e (1 - 2y) ef - 2i\ e y/yjl - y) l p e p /z u] 



(72) 



where {/iv} indicates symmetrization of indices and [fiv] indicates antisymmctrization. The fraction y is defined 
to be y = P2 ■ I/P2 ■ q ~ l~/q~, which in the lepton center of mass frame equals y = (1 ± cos# 2 )/2, where 9 2 
is the angle of hadron 2 with respect to the momentum of the incoming leptons. The contractions of specific 
tensor structures in the hadron tensor, given in Table 1, contain azimuthal angles inside the perpendicular 
plane defined with respect to Z^, defined to be the normalized perpendicular part of the lepton momentum Z, 



il/iQVyU^y))- 



l±- a± = — |aj_| cos 4 
e ± I/ /j_ M aj_ J y = |ajjsin^ a . 



(73) 
(74) 



V. INTEGRATION OVER TRANSVERSE PHOTON MOMENTUM 



In the next sections we will discuss explicit expressions for cross-sections. Instead of giving the complete 
cross-section, which can be obtained from the hadron tensor (Appendix B), we treat a number of special cases. 
In this section wc consider cross-sections integrated over all transverse momenta. 

After integration over the transverse momentum of the photon (or equivalently over the perpendic ular mo- 
mentum of hadron one Pi±_ = — ziq T ), the integrations over k T and p T in the hadron tensor (Eqs. (Bl) and 
(B2)) can be performed leading to 



f d 2 q T W^ = 12e 2 z lZ2 J2 e a 2 



M -9a 



D 1 D 1 -\ 1 \ 2 G 1 G 1 



(s&S& + g?S 1± .S a 



fli-ffi 
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TABLE I. Contractions of the lepton tensor L M „ with tensor structures appearing in the hadron tensor. 



-v 0- - y) \<*±\ \ b ±\ cos((j) a + 4> b ) 

y(l-y) \a±\ \b±\ sin(cp a + (/> b ) 
-(1 — 2i/)- v /y(l -y) \aj_\ cos(j> a 



til/ 

-STL 



a { /b» ± } ~(a ± -b ± )g» 






tjL u Xp 



(e^ax p b± a )g'T 



*<V? 



$ {m e ^>p , 



iaj^ 1 



(1 - 2y)^j/(l-y) |o±| sin<p a 

-A e (i -J/J |ttx| l&J-l Sin((j>b -<f>a) 



iz^e" ]p a 



Ae ^/y(l-2/)|a x | sin a 
Ae ^/y(l - y) [ax] cos<p a 



£ {/ x! 



Q 



A 2 
Sixp 



M 1 —G 1 +M 2 H 1 — 

z\ z 2 






Q 



«1 ^2 



Q 



: Ai 



M 2 Gi— + Mi— JTi 

^2 21 



Z fcj_ J2Xp 



Q 



D T # — 

M 2 D 1 —+M 1 — H 1 
z 2 z\ 



(75) 



and 



Jd 2 q T W^ = 12e 2 z lZ2 J2e a 2 



A1G1D1 - AaDiGi 



gl^c^] 



2z- 



IX 



£>T— #L 

M 1 — G 1 -M 2 H 1 — 
Z\ ' z 2 



+ 2i- 



AV"]Pq 

z e |_ oixp 



Q 



Mi— Di + MafTi— 

21 22 



2i- 



z o 2J _ 



Q 



-Ai 



M 2 C?i- 



2z- 



sMpo 
2 e ± o 2 xp 



M 2 D 



— - 

22 

Gt 

z 2 



Mx—Ex 

z\ 



Mi—Hx 

z\ 



(76) 



We have now included the summation over flavor indices and e a is the quark charge in units of e. The 
fragmentation functions are flavor dependent and only depend on the longitudinal momentum fractions, e.g. 
D\D\ = Df(zi) £)™(z 2 ). The result is expressed in terms of the fragmentation functions which survive 
the fc^-integration of the Dira c proj ections of the correlation functions (cf. Eqs. (E9f) to (J37J)): D\.Gi = 
Gil,Hi,E,El,Gt,H,Hl,Dt I35|j3q| . 

Note that the tilde functions arise naturally in the quark fragmentation region. The reason that this does 
not occur for the antiquark fragmentation is due to the non-symmetric choice of frame. This non-symmetric 
feature only shows up at subleading order. The leading order is symmetric, since ej_ acquires a minus sign, due 
to the interchange of the vectors n+ and n~. 

From the hadron tensors we easily arrive at the following expressions for the cross-sections, where we separate 
the cross-sections into parts for unpolarized (O) and polarized (L) leptons: 
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da°(e+e- -» h,h 2 X) 3a 2 ^ 2 J / x^rr 



dildzidz 



+ B(y) \S 1T \ \S 2T \ cos(<p Sl + <t>S 2 ) (ffi^i 

, x , v , , , x [2M 1 D T — 2M 2 H 
+ C(y)D(y) \S 1T \ sm(0 Sl ) [^^ D ^ + ~^ Hl J 2 

+ C(y)D(y) \S 2T \ sin^J (2^^ + ^^ 

+ C(»)D(y) A 2 |S 1T | cos^J ( 2 -^^G 1+ 2 -^ Hl ~i± 

\ Q z\ Q z 2 

- C(y)D{y) X, \S 2T \ cos(0 S2 ) (^G^ + ^^ ) } (77) 



and 



da L (e+e- -» WQ 3a 2 2 j C{y) ( , - 

a, a ^ x 



H-DO,) |S 2T | cos(<fe 2 ) f^D^ + ^1^ 
V Q ^2 Q zi 

, n ,„ , , , s /2Mi G T — 2M 2 I? 
+ £>(y)|5 1T | cos(0 Sl ) — i -^1 + -TT-ffi- 

- D(y) Ai S 2T smfe ~?r Gl n Hl 

\ Q z 2 Q Z\ 

+ D{y)X 2 \S 1T \ sin^J (W^-^Hi^V (78) 

where dSl = 2dyd(f> 1 , with (f> 1 giving the orientation of lj_, see Fig. pi Note that on the r.h.s. of the above 
equations the dependence on (f> 1 enters in the azimuthal angles, which are defined with respect to l^_, cf. Eqs. 
([73[) and (J74). We use the following factors: 

My) =\2^ v + y2 

B(y)=y(l-y), 
C(y) = 1 - 2p, 



D(y) = y/y(l-y). (79) 

The first three terms in Eq. (J77p coincide with the ones found in [B7J , if one neglects the contributions associated 
to Z exchange in their Eq. (45)7 One observes that besides these three leading contributions, one finds subleading 
single and double spin azimuthal asymmetries. 

To reduce the expression to the 1-particle inclusive cross-section, one must take the fragmentation functions 
for a quark fragmenting into a quark (see Appendix A) and sum over spins. Only Z?°(zi) survives and after 
summation over spins becomes a delta-function. We find for the one-hadron inclusive integrated cross-sections 
(using h as running index instead of 2 and realizing that D\ — D®): 

da°(e+e- -> hX) 3a 2 v-^ ?\ A , ,^ n , , „, s t>, s,„ , • ,, ,2M h D%{z h )\ 

dnd Zh l = W J2^\My)DU Zh ) + C(y)D(y)\S hT \ sm(fej ~^^\ (80) 

and 

da L (e+e- -+ hX) 3a 2 , v-^ ,f C(y), „.. ,„,,,„, ,, ,2M h GUz h )} 

-^mdz~ h " = W Xe 5 ea j " 2 XhG ^ + D ^ |5 -' cos ^ -cT^T 1 )- (81) 
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If the hadrons are unpolarized we find: 

da '(e + er — > hX) 3a 2 „. , v-^ 9 ^„, s ,„„% 

— 4nd£ — " = Q^ A ^ E e « 2 D ^)> ( 82 ) 

a, a 
a, a 

and d<r L = in both cases. Hence we find for the number of produced particles 

N h {z h ) = J2 e a 2 D<£{z h )/ J2 e a 2 , (84) 

N hlh2 ( Zl ,z 2 ) = Y, ea 2 D a l {z 1 )D'i{z 2 ) / £ e a 2 . (85) 

a,a a,a 

The case of Sh — gives the number of particles produced per spin degree of freedom, while the part proportional 
to Xh gives the contributions of produced hadrons with A^ = ±1. Thus the ratio of the part multiplying A^ and 
the Sh = result gives the longitudinal polarization of the produced hadrons, which must lie between —1 and 
+1. Similarly, the ratio of the part multiplying S hT and the Sh — result gives the transverse polarization, 
again a number lying between —1 and +1. In many cases the final state hadron will not be a stable particle, 
e.g. a A. In that case the final state (Nit for the case of a A) is used to determine the spin vector Sh pTJ. 

For the 1-particle inclusive cross-section we see one leading polarizing effect, namely for polarized leptons the 
longitudinal polarization of produced spin-1/2 particles is given by 

(longrtudmal polamataonj = -A e ^ ^T^I^T (86) 

At subleading order transverse polarization in the final state is induced given by 
transverse polarization \ _ . D(y) 2Mh 2-, a ,a e a t*T\ z h 



in lepton plane ) ^ Ae A(y) z h Q £ a , a ejD^Zh) ' (8?) 

transverse polarization \ C(y)D(y) 2Mh Z^ a ,s e °- D^(zh) 

(88) 



transverse to lepton plane J A{y) z h Q J2 a ,& e a 2 D1(z h ) ' 

where the lepton plane is spanned by I and P 2 . The in-plane polarization is proportional to the lepton polar- 
ization and is determined by the fragmentation function Gt ■ This function is the equivalent of the distribution 
function gx- An out-of-plane polarization is found for unpolarized leptons, determined by the time reversal odd 
fragmentation function Dt- The asymmetry ( p8| ) was first discussed by Lu p8| . 

For the 2-particle inclusive cross-section in which one hadron has spin 1/2, e.g. e + e~ — ► A7rAT, a longitudinal 
A polarization is induced, 

/, ,,, ,.,1 . C(y) E Q , a ea 2 Gr A (^r^2) 

(lon gl tudmal polanzat^n) = -A e ^ ^^a^^^ - ( 89 ) 

involving one polarized fragmentation function, namely G" - * . If both hadrons have spin-1/2 a correla- 
tion between the polarizations of the two hadrons exist. The correlated longitudinal polarization involves 
— A1A2 A(y) ^2 a a e a 2 GiGi; The correlated transverse polarization involves B(y) \S 1T \ \S 2T \ cos(<fisi + 
4>s 2 ) Eao 6 " 2 ^ 1 ^ 1 anc ^ provides a possibility to measure the transverse spin fragmentation function H\(z), 



the equivalent of the transverse spin distribution function h\ |37|. For the 2-particle inclusive cross-section there 
are several single spin asymmetries in unpolarized and polarized scattering, which among others give rise to 
twist 3 fragmentation functions Gt, Dt, Hl, H, E and El, in principle each with characteristic final state 
polarization, but suppressed by Mh/Q- 
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VI. LEADING ORDER ASYMMETRIES 

Instead of integrating out the q T -dependence, we will now focus on the fully differential cross-section, i.e., not 
integrated over transverse momentum (P\x = —Ziq T ). We will see that the transverse momentum dependent 
cross-sections contain asymmetries, which would vanish upon integration. Some of those asymmetries appear 
at leading order, to which we restrict in this section. The subleading results can be obtained from the hadron 
tensor in Appendix B in a similar way. 

First we consider the expression for both hadrons unpolarized: 



da°(e+e- -> h x h 2 X) _ 3a 2 2 2 



dSldz\dz 2 d 2 q T 



Q 2 



A{y)T[D l D l ]+B{y) cos(2&) T 



(2h-k T h 



''P r v — fC r p ' Pt 



MiM 2 



where we use the convolution notation 



T[DD] =J2e a 2 J d 2 k T d 2 p T 5 2 (p T + k T - q T )D a ( Zl , z 2 k^)D a (z 2 , z 2 p^), 



(90) 



(91) 



and da L = in this case. The angle <f>i is the azimuthal angle of h, see Fig. 1. So we find that the number of 
produced hadrons has an azimuthal dependence: 

N hl h 2 (zi,z 2 , q T , y) = z\z\ I T [DiDi] 



B(y) 
My) 



cos(2(^i) T 



Zi it ' IX, rp ft, ' Z/T~i Krp ' Z/T~i 



M X M 2 



E< 



(92) 



This asymmetry (the second term) has no analogue in the Drell-Yan process or semi-inclusive lepton-hadron 
scattering, since it involves a product of two time-reversal-odd functions. This new asymmetry goes with the 

same function H^ as appears in the Collins effect, multiplied with the similar time-reversal-odd function H 1 . 
We emphasize that this is a measurement in which no polarization of the produced hadrons is needed and the 
result is not suppressed by a factor of 1/Q. This in contrast to the cos(20) asymmetry found by Berger |2l| , 
which does not arise from time-reversal-odd functions. It is 1/Q 2 suppressed and also arises in other processes. 
Assuming for instance a Gaussian fc T -dependence of the functions, the convolutions can be evaluated. The 
number of produced hadrons would then be: 



N hlh2 (z u z 2 ,Q T , y) = G(Q T ; R) J^ ^ 2 I D a 1 {z l )D a 1 {z 2 ) 

a,a \, 



B(y) 
My) 



008(20!) 



QrpR 



M 1 M 2 R{R. 



lQ2 Ht a { Zl )Ht a {z 2 ) 



E^ 2 



(93) 



where R 2 = R\R]U{Ri + R$) an d Di(*i, k' T 2 ) = L>i(zi)i? 2 exp(-i? 2 fc 2 .)/ 7 rz 2 = Di(*i)0(|fc T |; Ri)/z 2 , etc. For 
details see Ref. Q. 

In case we consider the expression for hadron one polarized and hadron two unpolarized, we find the following 
additional terms: 



da°{e+e- -c h x h 2 X) 



ia 2 



2 2 

- 1 -9 



dfldzidz 2 d 2 q T Q 2 

-A(y)\S 1T \ sim>!-0 Sl )JF 

+B(y)\S 1T \ sin(30 x - <f> Sl ) F 



hk 



+ B(y)X 1 sin(2^i)^ 



I Z it' Krp ft' pi-p ~~~ rC n~< ' 'P r r 



Afi 



B(y)\S 1T \ anfa + fa) F 



hp T 



H 1L H 1 

MiM 2 



h x h\ 

Mi 



[4h-p T (h-k T ) 2 -2hk T k T p T - hp T fcf,J 2 M 2 M 



(94) 
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da L {e+er -> h x h 2 X) _ 3a 2 2 2 J x C(p) Ai ^ ^^ , 



dndz 1 dz 2 d 2 q T Q 2 \ 2 

Git-Di 



-A e ^|S 1T | 008(0! -0 Sl )^ 



/X ■ tZ^p ■ 



Mi 



(95) 



Again there is a term which has no analogue in semi-inclusive lepton-hadron scattering, namely the term with 

D^~ T , which can be seen by comparison with the result obtained in Rcf. p9| |. The term with H\H 1 is the analogue 
of the single- (transverse)-spin Collins effect. Note that it appears together with other single-transverse-spin 
asymmetries. 

Conversely, one can consider hadron two polarized and hadron one unpolarized, which may be simpler from 
the experimental point of view, because one does not need to measure the transverse momentum and the 
transverse polarization of the same hadron. We find similar expressions in which all single spin terms have, 
besides the obvious replacements, a sign change. 

The leading order double spin asymmetries can be found in Appendix C and are useful in for instance the 
case of e + e~ — ► AAX. 

We like to point out that the transverse momentum dependence of some of the functions can be directly probed 
in the situation where hadron two is taken to be a jet, which in this back-to-back jet situation is equivalent 
to analyzing the azimuthal structure of hadrons inside a jet. Only D\ remains and is a delta- function, so the 
convolutions can be evaluated exactly. In that case Eqs. (|90|), ( |94|) and (J9g) taken together yield: 



da(e+e -> h jet X) 3o , ^ - I . . 
d Sldz h <Pq T -~& Zh 2^ ea \ A[y > 



-A R « 



D1{z hl zlQ 2 T ) + \S hT \ sin(^ - d> Sl ) ^fDt£(z h ,z 2 h Q 2 T ) 



X h G a 1 {z hl zf l Q T ) + \S hT \ coa(4> h -<p Sl ) jfGlAzh, 4Qt) 

Mh 



(96) 



This result means that there is a transverse polarization transverse to the hadron plane, which is proportional 
to the function D^~ T , and a transverse polarization in the hadron plane, proportional to Git- 

So one sees that by measuring q T one can learn about the transverse momentum dependence of four functions 
in this particular case. There are no chiral-odd functions in this result, because they must be accompanied by 
a quark mass, which gives a result proportional to m/Q, so they are present in the subleading result. 

VII. WEIGHTED CROSS-SECTIONS 

The expressions in the previous section contain convolutions, which are not the objects of interest, rather we 
want the (universal) fragmentation functions depending on z and k T . At the end of the previous section we 
discussed a situation in which the transverse momentum dependence of some of the functions could be extracted 
from the analysis of one jet. Below we will outline a way to obtain instead of the full transverse momentum 
dependence, the fe|,-moments of the functions, defined as: 

F^(z 1 )=Jd 2 k' T ^yF(z 1 ,k' T % (97) 

for a generic fragmentation function F. The lowest moment is the familiar fc T -integrated fragmentation function. 
By constructing appropriately weighted cross-sections the convolutions result in products of such fey-moments. 
The same fe^-moments for instance show up in semi-inclusive lepton-hadron scattering, in that case multiplied 
by fey-moments of distribution functions |12j ] . In section 5 we have presented the hadron tensor and cross-section 
integrated over transverse photon momentum. A number of structures averaged out to zero, which are retained 
when the integration is weighted with an appropriate number of factors of q T . 

We find for the once-weighted cross-sections, where we again show only the leading results, for the case that 
hadron two is unpolarized: 



/, 2 / x d<r(e + e — ► hihiX) 2>a 2 ^— * 9 

d qT ( ^' a) dndz ld z 2 d 2 q T =-&2^ e ° |a 
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A(y) \S 1T \ sin(0 Sl -0 a ) [MxD^dA -B{y) \S lT \ sim> a + Sl ) [M 2 H 1 H 1 
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\S 1T \ cos{cf> Sl -cf> a ) UhG^D,] 1. 



(98) 



Constructing from this cross-section the weighted one-particle inclusive cross-section, by replacing D\ by a 
delta- function, and considering the specific case a = Zj_, one finds: 
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i-L(l) 



■-(1) 



These fc^-moments Z3 13 ; ; and G-j^ are related to the twist three functions Dt and Gt via 



£> T (z) = z s 



dz 



D#%) 



G T (z) = G 1 (z)~z 3 — 
dz 



G 1T (z) 



(100) 
(101) 



respectively [[12[ . An experimental verification of these relations by comparing the above cross-section to the 
one-particle inclusive results Eqs. (pfl) and (pTl) would be very interesting. 
The twice-weighted cross-section is: 
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(102) 



In particular, one can use (q T ■ li_) 2 , so one puts a = b = l± in the above equation (</> a = (f>b = 0), such that in 
case both hadrons are unpolarized 



d 2 q T fq T -/j_j 
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(103) 



Going back to the result for Nf H h 2 m Eq. (p2), one sees that weighting that result only with cos(2</>i) would 

produce a convolution of H^ and H 1 , while the result above shows that including appropriate factors of \q T \ 

produces a product of fc^-moments of fragmentation functions, in this case H 1 'H 1 . The fc^-moments can 
be used in other processes where they also occur. The above is an illustration of a general procedure. 

In Appendix D and E we give the integrated once and twice- weighted hadron tensors, respectively, in case 
both hadrons are polarized. The twice-weighted result is only given to leading order. 
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VIII. SUMMARY AND CONCLUSIONS 

We have presented the complete tree-level result up to order 1/Q for inclusive two-hadron production in 
electron-positron annihilation. We consider the situation where the two hadrons belong to different, back-to- 
back jets. Polarization in the initial and final states is included for the case of spin-1/2 hadrons. In case of 
spinless hadrons one will focus on the ones that are produced most abundantly, like 7r's and -ftT's, which also 
serve to study flavor dependence of fragmentation functions (see for instance |4(|). For the case of spin-1/2 
hadrons A's seem most appropriate due to their self-analyzing decays (see for instance 37 411] )• Hadrons with 



higher spin, like p's, are not considered, because in that case a spin vector is not sufficient to describe the spin 
states. 

We have restricted ourselves to the case of one photon exchange, since we are interested in power corrections 
which are, most likely, negligible in regions of Q 2 where the annihilation into Z bosons becomes important, 
i.e., LEP energies. In forthcoming work we will investigate the inclusion of Z's in the leading part of our result 
taking into account transverse momentum. 

We have worked in a diagrammatic approach based upon analogy to the one developed by Ellis, Furmanski 
and Petronzio for DIS. In this approach soft parts of a scattering process are treated as hadronic matrix elements 
of non-local operators. The soft parts occuring in the process under consideration are given by quark-quark and 
quark-gluon-quark correlation functions. The latter are necessary to achieve electromagnetic gauge invariance 
and are related to the quark-quark ones by use of the equations of motion. All soft non-perturbative physics 
is then parametrized by a set of fragmentation functions. The leading ones can be interpreted as quark decay 
functions. 

We have done our calculations at tree-level under the assumption that collinear divergences could be absorbed 
in the fragmentation functions and that factorization consequently holds. Moreover, for the interpretation of 
our results one should keep in mind that loop corrections, in general, will lead to finite order (a s ) "-corrections 
effecting the magnitude of observables and may also lead to non-zero contributions to observables not present 
at leading order (a well-known example is Wl in DIS). 

Our results include among others the following: 

• We have considered the cross-sections integrated over transverse momenta of the produced hadrons for 
both, polarized and unpolarized beams up to subleading order. Our result contains a few terms, which 
have been found previously in a leading order analysis for unpolarized e + e~ annihilation p7| . 

• In particular, we have focussed on the information obtainable by observing transverse momentum of one 
of the produced hadrons (defined either relative to a jet-axis or relative to the momentum of a hadron 
in the second jet). Although cross-sections differential in transverse momentum are not easy measurable, 
they are of particular interest, since they contain leading order asymmetries, which would vanish upon 
integration. 

We have found a number of new unpolarized, single and double spin asymmetries. Often they have 
no analogues in (semi-inclusive) DIS or the Drcll-Yan process, since they involve products - or to be 
more specific, convolutions - of two time-reversal odd fragmentation functions. In particular, the cos(20) 
dependence discussed in Sect. 6 is most likely measurable, since it is not suppressed by powers of 1/Q and 
does not involve polarization, neither of the beams nor of the final states. 

One-hadron inclusive measurements supplied with the additional determination of the jet axis gives direct 
access to the transverse momentum dependence of some of the fragmentation functions. 

• We have discussed how convolutions of fragmentation functions can be converted into products of their 
fcy-moments. This is achieved by appropriate weighting the integration over the transverse momentum 
dependence, in the spirit of a Fourier analysis. This is another way of retaining asymmetries, which would 
vanish upon (non-weighted) integration. 

As a final note, it can be seen from our results which extra asymmetries will show up in, for instance, the 
Drcll-Yan process if one allows for time-reversal odd distribution functions. Such a single spin asymmetry 
is discussed in p2| . However, the presence of time-reversal odd distribution functions would require some 
factorization breaking mechanism, like the one discussed in p3| . 

In conclusion, we emphasize the possibilities provided by measuring azimuthal asymmetries and higher twist 
annihilation in order to learn more about the structure of hadrons. 
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APPENDIX A: THE FRAGMENTATION FUNCTIONS FOR A QUARK FRAGMENTING INTO A 

QUARK 

We consider the correlation function for the case of a quark (with momentum k) fragmenting into a quark 
(with momentum p and spin s), given by 



Sij(p,s;k) = u l (k,s)u j (k,s)6 i (k - p) = -((# + m)(l + 75^}) ij S 4 (k - p), 



(Al) 



where the momentum and spin of the quark are parametrized as 
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-, k T 



X q k 



mX q 
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k T ■ s qT 



y\n Krp 



2mk-' qT 






(A2) 
(A3) 



in terms of a quark lightcone helicity X q and a quark lightconc transverse polarization s T . The projections 
become for twist two 



5 b-] {k )= l -5{z-l)d\k T -p T ), 

5^~^\k)= l -X q S(z-l)5 2 {k T -p T ), 
S^-^(k) = lsi T S(z-l)6 2 (k T -p T ), 



(A4) 
(A5) 
(A6) 



where z = p jk . For twist three we get 



S [1] (k) = ^-S(z~l)S 2 (k T -p T ), 

si-n ik) = 3-5(z-i)s 2 (k T -p T ), 
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APPENDIX B: THE COMPLETE EXPRESSION FOR THE HADRON TENSOR 



The full expressions for the symmetric and antisymmetric parts of the hadron tensor are (expressed in the 
perpendicular frame defined in section 2) 



Wg" = 12e 2 ziz 2 / d 2 k T d 2 p T 5 2 {p T + k T - q T )< 



-9? 



D\D\ — Gi s Gi s 



e^k±pSi±o 

Ml 



D± T Di 



pa c 

C± P±p0 2 ±a „ —± 



M 2 



-DxD 



IT 



20 



k±p± Si±-S 2 ± - p±-S i± fc±-5 2 ± n± j^± ' 

L) 1T D 1T 



Ap. qv} , av 



— [S 1A _S 2A _ + g ± Su_-S 2 ± ) H\tH\t — 



M 1 M 2 

— k^p2 +gT k ^-PA 



MiM 2 



H ls H ls + H l H l 



M H ^ HlT M 2 HlTHls 

^-^- \H U H X H, H u 

k[^S 2±p + S { 2 ^k ± pM P Si± p + S&<Q>p± p -l 
^1 -Hit H T777 ti VT tl 1 



2Mi 



2M 2 



rH'? 



Z^- S 1JL -S 



-£>i - 



•^2J 



Z\ 



M 2 DJ T ^-\ 2 P± - Sl± 

M x 1T z 2 M 1 



D 



IT" 



z 2 



, fcj.-p± gi±-Sg± - P±-S 1± k±-S 21 _ . — ± 



M 1 M 2 



Zl 



H rp — — 



M 2 



M 2 21 «1 Mi l Z2 z 2 



#s 



# 



Z Pi 



Q 



z 2 



ir k\ s 1± -s 2± 



2MiM 2 



Dyr^iT "I" -V 



fc ± -S 2± Df ■ 



M, 



Dt^I- 



-D 1T + — 5i_L • S 2 ± —Dit 
M 2 Z\ M 2 Z\ 



z 2 Mi z 2 z 2 M 2 \zi zi 



z D 1± 



+ A ; 



^■^-D^ + ^k^.S^D^-M^Gu 



Mi 



z 2 Mi 



— — H 1 — M 2 Hit k±-S 2 ± Hit 

M 2 z\ 



He 

z 2 



_|_ Dt . , G'rr — 

z 2 z\ 

Z\ 



&S"J - 



Ai 



k±p± Dj; — i , P±-Si±k\ :-i_ Mi 



M 2 zi 



-D 



D T -i 



IT 



2M!M 2 ^ 1T ^ 1T ~ M- 2 P± - Sl± ~ DlT 



M 2 G ls ^ + ^Ht^ 

Z2 Ml Z2 



p_L • SiJ. i7 it — + Mi —H 1T 

z 2 z\ 



r { V; }p fc 



-i-P 



X^.E^D 



i._ ,-, 1 -D fe-L-Sij. ,, 



-21 



Mi 



IT" 



~2 



Mi 



m ,— p±-Sii_ H T —j_ 
— Jrl 1 Lr ls — ti x 

Mi M 2 z\ 



Df T D 1 

Mi \ ls z 2 1 z 2 



z^e ± rr p ±p 



A 2 Z?i 



D L k ± -S 21 _ D^-± k^-S 21 _ ,H T Mi H 



--■2 



M 2 zi 



D 1T - 



M, 



H 



1 z 2 M 2 \ Zl 



Z\ 



z,{V c "}Pq 

Z 6j_ 0\± p 



Dt-^t k±p± Rk — ± fei-pi „ 1 D k 2 , ^ — i? 

M 1 —D 1 + y —^H, + ^ Di T + ^^Di T D x + M 2 H 1T — 

zi M 2 zi x Mi 1T z 2 2Mi 1T 1T 2 2 



z e_|_ o 2 ± p 



_ Dt if — k±pi D — ± fei -pi 1 _ff T 
M 2 D 1 —!- + M X —E XT + ,, D 1T + ^ -ffi — 

z 2 «i M 2 zi 1T Mi x z 2 



(Bl) 



21 



and 
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APPENDIX C: DOUBLE SPIN ASYMMETRIES 



In this Appendix we give the azimuthal dependences of double spin asymmetries, as can be observed, for in- 
stance, in AA production by determination of the polarizations of both observed hadrons. The spin-independent 
and single-spin dependent parts of the cross-section are given in Eqs. (|90|), (p4) and (95). 
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APPENDIX D: INTEGRATED ONCE- WEIGHTED HADRON TENSOR 

We display the hadron tensor weighted with the factor ( q T ■ a ) and integrated over the transverse photon 
momentum. The vector a is an arbitrary vector like, for instance, l±. 
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APPENDIX E: INTEGRATED TWICE- WEIGHTED HADRON TENSOR 



We display only the leading terms of the hadron tensor weighted with two factors (q T -a) (q T -b) and 
integrated over q T . 
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